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We tested the Horava Lifshitz (HL) quantum gravity model by using the Lii, Mei and Pope 
solutions for primordial black holes (PBHs) and the observational upper bounds of the PBH density 
parameters. We found that, although the HL model is severely constrained, it is not ruled out. 
When our analysis is combined with that of Dutta and Saridakis the observed value of the density 
parameter Upbh might rise by several percent as the running energy parameter A increases. 
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I. INTRODUCTION 

Horava L lj has applied the anisotropic scaling that Lif- 
shitz developed^ for quantum gravity and established 
the power-counting renormalizability of the theory in 
3+1 dimensional space-time. This model is usually called 
Horava Lifshitz (HL) gravity. Cosmological solutions 
of HL gravity for the Friedmann-Lemaitre-Robertson- 
Walker metric have been described Lii, Mei and Pope 
(LMP) ||. Although there has also been a lot of other 
work on HL gravity Q, the LMP solution is the most 
basic one, and therefore our discussion is based around 
it. The LMP solutions are based on conditions of relaxed 
projectabiiity and of detailed balance. Recently, Gong et 
al.[5j have shown using linear cosmological perturbation 
theory with projectabiiity that HL gravity is reduced by 
the Faddeev and Jackiw condition[6j and that there is a 
scalar ghost mode for i < A < 1, where A is the running 
energy parameter in the HL model. 

In 1971, Hawking showed that black holes should 
have been generated with various masses in the early 
universe [7|- Nowadays, such black holes are called pri- 
mordial black holes (PBHs), and their masses M are esti- 
mated from the Hubble equation as M ~ 10 15 ( 10 ^ 23 j g, 

where tf is the current time from the big bang. The black 
holes emit energy according to the Stefan-Boltzmann law 
and evaporate. Because the time necessary for evap- 
oration is t ev ~ 10 10 (j^rs-g) 3 years, black holes that 
are lighter than 10 15 <? have completely evaporated via 
Hawking radiation before the present t/(— 10 10 yr)[l| 
and only heavier PBHs remain. In Q, the observational 
values are given between 35 MeV and 175 MeV for 7- 
rays, Vt PBH < (7.6 ± 2.6) x lfr 9 ^ 2 ~ 1.4 x 1CT 8 , 
with the Hubble expansion rate h = 0.73^q'o3 ~ 0.73, 
and between 30 MeV and 120 GeV for 7-rays, SIpbh < 
(5.1 ± 1.3) x 10~ 9 h- 2 ~ 9.5 x 10~ 9 . The observational 
upper limit on the PBH density parameter from mea- 
surements of Hawking radiation is @ 

VPBH < 10~ 8 . (1) 

However, PBHs are as yet unidentified stellar objects [To| 
and may constitute all the dark matter [llj]. Regrettably, 



the nature of dark matter and dark energy is not yet un- 
derstood, but here we assume that the density parameter 
of the PBH fl p bh contains the density parameters of the 
dark energy VLee or the dark matter Udm- The values 
of the various density parameters from WMAP [T^] are 
n b ~ 0.0449 ± 0.0028 for baryons, fl DM ^ 0.222 ± 0.026 
for the dark matter and £Ide — 0.734±0.029 for the dark 
energy. 

In HL gravity the (running) energy parameter A lies in 
the range | < A < 1. The present value of the energy 
parameter is | Ao — 1 |< 0.002, based on observations 
of supernova, the baryon acoustic oscillation (BAO) and 
the cosmic microwave background (CMB) [1J]. However, 
because we are only discussing PBHs and the focus of our 
research is different from theirs, such a tight bound on 
Ao may not be appropriate. Note also that many works 
have been carried out that focus on dark energy in the 
HL modelpj]. 

This rest of this article is arranged as follows. In the 
next section we review HL gravity and the LMP solu- 
tions. In the following section, we compare the LMP 
solutions with current observational PBH data. Finally, 
we present our conclusions. The values in SI units of 
some cosmological constants are given in an appendix. 



II. HORAVA LIFSHITZ AND LMP SOLUTIONS 

HL gravity [l| has the anisotropic scaling property x — > 
bx, t — ¥ b z t, with critical exponent z. The theory is 
Lorentz invariant in the long distance, IR limit (z = 1) 
but it is not Lorentz invariant in the short distance, UV 
limit (z = 3). 

The Arnowitt-Deser-Misner (ADM) decomposition of 
4-dimensional space-time is given by 



ds 2 



c 2 N 2 dt 2 + g tj (dx 1 - cN l dt) (dx j - cN j dt) ,(2) 



where N is the lapse field, N 1 is the shift field and gij 
is the spatial metric. The scale dimensions of physical 
quantities in units of mass are [Gn] — —2, [N] — [gij] — 
0, and [Ni] = z — 1. For simplicity, we assume that N 
depends only on time t (i.e., projectabiiity). The ADM 
decomposition of Einstein-Hilbert (EH) action is 
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where Gn is Newton's gravitational constant and 
Kij is the extrinsic curvature defined by — 
W (dij ~ ^iNj — VjNi). The original action of 
Hofava[l] is 



Shl = J dtd?x{L +Ci), 

Co = yfgN |^ {KijKV - XK 2 
n 2 p 2 (A W R-3A W 2 ) \ 



(4) 



8(1 - 3A) 

acV(1-4A) 
32(1 - 3A) 
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(5) 
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where X,K,fi,w and Ay^ are constant parameters, and 

Cy is the Cotton tensor C 13 = e m V k (r{ - \R5fY 

Note that only w controls the strength of the interac- 
tion and the free state corresponds to the limit w — > 0. 
Also note that we can identify A — > 1 as the IR limit and 
A — » 1/3 as the UV limit. By comparing Q with © the 
following relations are obtained. 
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A = -Aw, 



(7) 



where A is the cosmological constant. In 3+1 dimensions, 
the correspondence of the HL action to the Einstein- 
Hilbert (EH) action demands A = 1 (i.e., the IR limit). 

If we demand that the speed of light is positive in ||7J) , 
then Ayy is negative. The analytic continuation fi —> 
ip and w 2 — > —iw 2 yields positive Aw- Then, we can 

describe the speed of light asc= f° r ^ > §■ 

The Friedmann-Lemaitre- Robertson- Walker form is 



ds 2 



c 2 dt 2 +a 2 (t) 



dr 2 



1 - Kr 2 



+ sin 2 6># 2 ) . (8) 



where K = — 1,0, 1 corresponds to a closed, flat or open 
universe, respectively, and a(t) is the scale factor. For 
this metric the Friedmann equation (in units where c = 
1) is 



3A- 1 

Aw (8ttGn)p 
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, (9) 



where p = p m + p r is the sum of the energy density of 
matter (p m ) and the energy density of radiation (p r ). 
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where the energy density of matter p m consists of the 
baryonic energy density pi, and the dark matter energy 
density pdm, the total energy density of radiation p r 
consists of the photon energy density p 1 and the muonic 
energy density p^ (i.e., p r = p 1 + p^ = 1.69p 7 [15] ), and 
H = ~ is the Hubble parameter. On the right hand side 
of (JTUJ), we identify the first term as the matter density 
parameter fi m , the second term as the radiation density 
parameter fi r and the third term as the curvature density 
parameter fix- Consequently we obtain 



K 



(3A- l)H 2 



K 2 



Awa 4 



A 



w 



(11) 



The right hand side of this equation is called the dark ra- 
diation [TJ|, [I?} ■ The dark radiation contributes negative 
energy density and the greater part of this contribution 
disappears if the curvature K is zero. 

We make the following assumptions. The density pa- 
rameter of PBHs takes its maximum value QpBH,max if it 
includes the density parameters of the dark matter fi^M 
and dark energy Qde- So 

QpBH,max < 0.96 = fi_D_E + fi/\ + &DM + fir- (12) 

From current observations (via Hawking radiation) [l8j 
(also see 0, [Tj| ) , it is thought that the present density 
parameter for PBHs Qpbh satisfies Qpbh < 10~ 8 (HI- 
Because {1} is based on data from radiation observations, 
SlpMo m ay contain the dark radiation term (i.e., CIde)- 
Thus, we obtain 
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(14) 



(3A - l)H 2 \A Wo a 4 

Because Aw is small, we can neglect second terms of 
(fTB")) and (fT4")) . Then, assuming K ~ K, Aw — Aw and 

a = 1, 

c 2 K 2 



fi 
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When we take the ratio of (fT5")) and (fTl?)) . we obtain 
fipsH 3Aq — 1 HqC/ 



fi 



PBH a 



3A - 1 H 2 a 4 cl 



0.96 x 10 5 



10 s 



(15) 
(16) 

(17) 



III. COMPARISON WITH PBH DATA 

Reintroducing the speed of light explicitly, the den- 
sity parameter is defined by the Friedmann equation @ 



Substituting ( ± J = ^Aof^ i n t JTTJ) gives 



Hq ( 3Aq — 1 



H 2 a 4 V 3A - 1 



10 8 . 



(18) 
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Consequently, 



That is, 



H 2 ~10 8 H 2 a 4 ' '' !A ' 



= 10 8 dV 



3A -1 
3A- 1 
3A -1 



Introducing r 2 = 10 8 , we have 



Hn 



„2 - 2 2 



( 3A- 1 
V3A -1 
From ([9]) the scale factor is 
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(3A-l)a 2 (3\-l)A w a 4: 
where p = p m + p r . Now, we define 

A w 2(8irG N )p 



A = 
B = 
C = 



3A-1 3(3A-1)' 
2K 



(3A-1)' 
(3A - 1)A W 



(21) 

(22) 
(23) 
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with the present values being denoted by the subscript 
0. Then (JUJ) becomes 



a = Aa — B + C- 



(25) 



which can be solved easily. This differential equation has 
two types of solution for the present Hubble parameter 
H : 



(I) 
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(27) 



where ao and /3o are the constants of integration. With 
assumption that the extrinsic curvature K is zero, we 
have Bo — Cq = 0. Then the constants of integration are 
given by 



a \K=o = ±2e-^ to ^Ao~, 
/3 | K =o = ±2e v ^ to v /^ 
and the present Hubble parameters are 

h-Wi l~T~ ( 3A - 1 



V3A -1 
3A- 1 



3An-l 



(28) 
(29) 

(30) 
(31) 
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PO = Pm + P la - Pm = 0.24 X 
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(32) 



(33) 



Then H 2 becomes 



2 _ ^(3A-1) 2 r 2 



(3A - 1) 

-36H /„2 



{A^-fO^tf, 2 }. (34) 



With = 7.77 x 10 - 36 [l/s 2 ] and H Q = 4 . 35x l i7 [a] (see 
Appendix), we obtain 



(3A - l) 3 A 



w 



r 2 (3A-l) 2 H 2 



0.48 = 1.95. (35) 



Thus, 



A = 



1 , 1 /(3A -1) 3 



3 3V 1.95 x r 2 



(36) 



In the HL theory the energy parameter A is in the range 



| < A < 1, so we obtain 



1 1 + ^7.80 x r 2 
3 < Ao < - . 



(37) 



The physical interpretation of ([57)1 is discussed in the 
next section. 



IV. CONCLUSIONS 

If 0.96 is substituted for ilpBH.max and 10 -8 is substi- 
tuted for QpBHo: (|37)) becomes 



1 



< A < 307. 



(38) 



Although this range contains the parameter range < 
Ao < 1) of the HL model, it is clear that HL gravity is 
considerably more restrictive. 

Inversely, when the limitations of the HL model are 
imposed on (|37[) we obtain 



PBH a — QpBH 



(39) 



This is consistent with the result of Frampton et al. [Tl| . 

If the black holes do not contribute to the density pa- 
rameter (s), the present energy parameter Ao in (|37l) be- 
comes 1/3 (i.e., we attain the UV limit). In the HL the- 
ory the general theory of relativity (GR) is not formed 
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in the UV limit, but is recovered in the IR limit. Thus, 
A — 1/3 is inconsistent with the status of GR of today. 

It is interesting to combine our analysis with that of 
Dutta and Saridakis 13]. If (|37|) is combined with their 
result | Aq — 1| < 0.002, it becomes 



current observational upper bound £Ipbh < 10~ 8 and 
the assumption QpBH,max — ^de/dm are correct, then 
the HL model is doubtful. (3) If HL gravity and the rela- 
tion (|39[) are correct, then the observational flpBH m ay 
increase slightly as A increases. 



1.02 X Q P BH < ttpBH,max < 1-03 X QpBH - (40) 

This shows that QpBH,max can be larger than present 
observed value flpBH by few percent. This is consistent 
with the prediction in the HL theory that the horizon of 
a black hole appears as A increases [H 0. A PBH that 
has not been observed up to now will be observed in the 
future because of an increase in A. 

We thus arrive at the following three possible conclu- 
sions. (1) If the HL theory and the current observational 
upper bound flpBH < 10~ 8 are correct, then the as- 



Appendix: Cosmological constants 



The present Hubble parameter is observed to be 



H 
71.0 

M. 



71.0 ± 2.5[km/(sM pc )} [12]. So we take H 



kin 



4.35X10"M ' Where 



SH^j = 2 - 30 >< 10-1 
,,, - 3.09 x 10 22 [m]. 
Generally the cosmological constant, i.e., the density 
parameter of the vacuum (dark) energy, is given by 



\ 



= 0.734. So we obtain A = 3 x 0.734 x = 



sumption flpBH,i 



fl 



DEI DM 



is doubtful. (2) If the 1.16 x 1CT 35 [i] and A\ 



3 J1 



7.77 x 10 



-36 
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